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Abstract. Light forces induced by scattering and absorption in elastic dielectrics lead 
Oh 1 to local density modulations and deformations. These perturbations in turn modify 

light propagation in the medium and generate an intricate nonlinear response. We 
generalise an analytic approach where light propagation in one-dimensional media 
of inhomogeneous density is modelled as a result of multiple scattering between 
J>-^ polarizable slices. Using the Maxwell stress tensor formalism we compute the local 

optical forces and iteratively approach self-consistent density distributions where the 
elastic back-action balances gradient- and scattering forces. For an optically trapped 
dielectric we derive the nonlinear dependence of trap position, stiffness and total 
■ deformation on the object's size and field configuration. Generally trapping is enhanced 

, by deformation, which exhibits a periodic change between stretching and compression, 

t — ' This strongly deviates from qualitative expectations based on the change of photon 

as : momentum of light crossing the surface of a dielectric. We conclude that optical forces 

have to be treated as volumetric forces and that a description using the change of 
photon momentum at the surface of a medium is inappropriate. 
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v~j . 1. Introduction 

rS 

b 

As light carries momentum besides energy, its propagation through a polarizable medium 
is accompanied by forces. Although the momentum of a single light quantum is very 
small, laser light can generate appreciable forces on the microscopic scale. Optical 
forces are nowadays routinely used to manipulate and trap particles ranging from single 
atoms and molecules [U El |3] to plastic beads, biological cells or microbes up to the 
size of tens of micrometres [31 El E] . The mechanical motion of even larger objects 
such as silica mircodisks or suspended mirrors has been damped and cooled by light 
forces E] . While most of the existing work targets the overall effect on the centre of 
mass of the particles, it has been shown by us as well as by other groups that these forces 
do not act homogeneously but exhibit distinct patterns within the medium [HI [101 E] ■ 
For any elastic medium this leads to local compression or stretching. Of course the 
modified density also changes the local refractive index and light propagation, which 
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Figure 1. Schematic illustration of the interaction between optical forces and local 
deformations within elastic media. 



again alters the forces as displayed schematically in figure [TJ The resulting coupled 
complex evolution thus obviously requires self-consistent models and solutions [12j. In 
addition, as the light mediated interaction is inherently long range, even a small but 
periodic variation of the refractive index can have a very large overall collective effect 
coupling distant areas over a large volume. 

This work is organised as follows: In section [2] we first present the basic scattering 
approach to treat the light propagation in an inhomogeneous refractive medium and 
use a previously developed formalism based on the free space Maxwell stress tensor to 
calculate the corresponding local force distribution (section [3]). This method is then 
used in section H] to develop an iterative scheme to calculate the steady state density 
and field distribution as a function of geometry and field intensity. In section [5] we 
discuss essential physical consequences predicted by our model at the hand of numerical 
examples. Finally in section [6] these results are set against common calculations of the 
total deformation at hand of the change of photon momentum at an interface between 
to dielectrics. 



2. Multiple scattering model of light propagation in inhomogeneous media 

The effective light propagation in a medium can be seen as the result of multiple 
individual scattering processes, which in general requires intricate numerical treatments, 
if one cannot make use of material symmetries. Here we restrict ourselves to the simple 
but still nontrivial case of two incoming counterpropagating plane waves in a transversely 
homogeneous and linearly polarizable medium. In this limit only forward and backward 
scattering add up phase coherently, while all amplitudes for transverse scattering average 
out. From the viewpoint of the forward and backward propagation directions, transverse 
scattering thus can just be added to an effective absorption rate in the medium. This 
is certainly not perfectly fulfilled in an actual setup, but still can be expected to give 
the correct qualitative behaviour, as long as the transverse extensions are much larger 
than the wavelength of the light. A more realistic treatment, e.g. in terms of Gaussian 
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transverse beams, is possible, but greatly complicates the model and would obscure 
many interesting physical phenomena found in this simple approach. 

Restricting the dynamics to the forward and backward scattering amplitudes along 
the propagation directions gives a simple and tractable model for our medium via a one 
dimensional array of N thin slices at positions x\, . . . ,xn- Here the spatial behaviour 
of the electric field E(ar, t) = Re[E(x) exp(— iut)]e y is determined by a ID Helmholtz 
equation p31 Q31 EE] 

TV 

(d 2 x + k 2 )E(x) = -2k(E(x) S (x ~ Xj). (1) 

i=i 

The field-induced polarisation density at each slice then is P(x) = anxE(x) J2f=i $(. x ~ 
Xj), where we introduced the dimensionless coupling parameter ( = krjAOt/ '(2sq) 
proportional to the atomic polarizability a and the areal particle density rj^ within 
the slice. Eq is the vacuum permittivity and k = to/c the wave number of the optical 
field. Note that we assumed here that the dipoles in each slice can simply be added up 
coherently for scattering along the propagation direction. As illustrated in figure [2, the 
equation above is satisfied by interconnected plane wave solutions [13], here denoted as 

Ej(x) := C,-e ifc(x J '- ] + Pjv :/ ' : '' = A J+1 e ik{x -^ +l) + B J+1 e- ik{x -^ +l) , (2) 

for Xj < x < Xj+x- The amplitudes left and right of a material slice (beam splitter) at 
position connected via 



Dlj V _ iC ^ )KBjl (3) 

The amplitudes (Aj,Bj) and (C 3 -_i, Dj-i) are coupled by a simple propagation matrix, 
i.e. (Aj,Bj) T = Pdj (Cj-i, Dj_i) T with P d . := diag(exp(iA;c?j), exp(— ikdj)), with the 
distance dj := Xj+i — Xj, j = 1, . . . , N — 1. 

Therefore, the amplitudes to the left of the (j + l) th slice are obtained by a simple 
multiplication of the previous transfer matrices, 

(s'+J = P ^ M BS---Prf 2 MBsPd 1 M B s (Jfj ■ (4) 

The amplitudes Ax and are determined by the amplitudes and phases of waves 
coming in from the left (i. e. — oo) and from the right (+oo), respectively, and constitute 
boundary conditions on the solutions for the Helmholtz equation (CQ). B\ and Cn are 
obtained by computing the total reflection and transmission amplitudes via 

^( t2 ^ irr ^ =M B sP rfjv _ 1 MBs---Pd 2 M B sP (il MBs, 

Bx = nAx + tD N and C N = tA x + r r D N . 

Note that the reflection coefficients for left or right incidence on an inhomogeneous setup 
usually do not coincide, i. e. r\ ^ r r , but the transmission amplitude t is independent of 
the direction of propagation. More details on the properties of these generalised transfer 



matrices are given in Appendix A 




Figure 2. Schematic illustration of the displacement and deformation process i4i = 
x+u(x). The initial medium (left figure) occupies the space [0, L], the slices are marked 
with dotted lines separated by do = L/{N — 1). In this case, the field generated by 
multiple scattering by the beam splitters (blue curve, cf. ([2])) reproduces the solution 
for a homogeneous medium with refractive index n (red curve) , if the coupling is chosen 
as in ([6]) . On the right hand side we see the displaced medium with irregularly spaced 
slices of the same coupling £ and the resulting electric field. The background shading 
illustrates the change in the distances, i. e. the strain u'(x) = — (p(x) — p)/p, with dark 
colours indicating regions of higher density. 



For equally spaced, thin polarizable slices we set Xj = (j — l)d , such that x = 
and xn = (N — l)d =: L and (J4j) simplifies to (Aj+i, Bj+i) T = (Ai, Bi) T , with 
Th := P^Mbs- In an earlier work [TT] we showed that choosing a uniform distance do 
between the slices and setting the coupling parameter 

cos(kd Q ) — cos(nkd ) 
sm(kd ) 

leads to the same optical fields as found inside a medium with refractive index n. A 
sufficiently dense array of beam splitters with spacing do = L/(N — 1) is then in the 
limit iV — > oo indistinguishable from a homogeneous medium of refractive index n and 
length L. 

A decisive step in this work, which allows us to account for local material density 
variations, is the introduction of a local displacement variable u(x) 

Xj i — y xj = Xj + u(xj), j = 1, . . . , N. (7) 

As illustrated schematically in figure |2l such shifts alter the local fields as well as the 
total reflection and transmission properties of the object. 
The distances between the slices then change as 

dj — d = Xj+\ — Xj — d = u(xj + i) — u(xj) =: Aj. (8) 
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A continuous limit can be consistently denned via u(xj) — > u(x) for x 6 [0, L] to obtain 
lira ^ = lim ^ + ^ = ^ = «'(*). (9) 

In analogy with the theory of elastic deformations, we call u' strain or deformation [THJ 
US], and the relative change in the initially homogeneous local material density p simply 
reads 

~ P{X) ~ P = -u\x). (10) 
P 

Let us here comment on the notation we will use for the rest of this work. As defined 
in the paragraph above (jSJ), our coordinates are chosen such that the unperturbed 
medium occupies the region [0, L\. Introducing a displacement u then shifts the object 
to [u(0), L + u(L)], with L + u(L) — u(0) =: L. But to ease notation, all the quantities 
such as the electric field strength E or force F shall remain defined with respect to the 
original position such that e. g. E(0) [E(L)] always marks the field at the left [right] 
edge of the medium. The amplitudes at the boundaries then have to be adjusted with 
corresponding phases, cf. (122|) . This, however, is relevant for mathematical formulations 
only, physical discussions and figures are unaffected by this detail. In © we introduced 
a tilde to distinguish the shifted Xj from the original Xj. For most other quantities such 
as the fields or forces, we will omit this tedious notation. Only the changed length L, the 
inhomogeneous density p ( |T0|) and refractive index n ( |T9l still have to be distinguished 
from their original values L, p and n. 

As mentioned before, defining the coupling ( as in (Q ensures that the solutions of 
the wave equation (CQ) agree with the field inside a homogeneous dielectric at positions 
Xj = (j — l)d , if the fields are assumed to agree at x\. In the continuous limit iV — )■ oo, 
the latter requirement is always fulfilled Interestingly we still preserve this feature 
for a model with displaced slices, if we choose the following approach: 

Let, as in ((2]), Ej(x) denote the plain wave solution of the Helmholtz equation (CQ) 

and 

E„ (.r) (;^'^ ,K r 1 ^ ■ Iljv 'V for Xj < x < x j+1 (11) 

denote a field defined in the same region, but with a refractive index rij. To obtain the 
desired equivalence between a stratified dielectric and a set of irregularly spaced slices, 
we assume for any given j e {1, . . . , iV — 2} 

lim Ej (x) = lim E n . (x) (12) 
and demand that with E{^y) = \im x ^ y E(x), 

E j(txj+i) = E j+1 {ix j+1 ), E'jtfxj+t) = E' j+1 (lx j+1 ) + 2k(E j+1 (lx j+ i), (13) 

E n] (t^'+i) = E nj+1 (ix j+1 ), E' n . (tx j+1 ) = E' nj+i (\x j+1 ), (14) 

Ej(txj+i) = E nj tfx j+1 ), E j+1 tfx j+2 ) = E n . +1 tfx j+2 ). (15) 

The first line shows the conditions that Ej and Ej + i are solutions of the Helmholtz 
equation ([1]), cf. ([3]) or [13J, the second line denotes Fresnel's equations for the transition 
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between two dielectrics. In the third line, finally, we demand that the plane wave 
solutions of §B) agree with the fields inside the dielectrics at positions Xj + i and Xj +2 - 
This leads to the required, successive coupling between £, the distances dj = Xj — Xj-\ 
and dj + i, and some indices rij, rij + x. 

Solving (Tl3l) -( fT5l) under the assumption (fl2l) and demanding solutions independent 
of the field amplitudes results in two conditions, for j — 1, . . . , N — 1 
rij sin(kdj) n,j+\ sin(kdj + i) 



sin(knjdj) sin(knj + id 



cos(kdj) — cos(rijkdj) cos(kdj + i) — cos(rij + ikd 



+ 



(16) 
(17) 

sm(kdj) sm(kdj + i) J 

One can easily check that these conditions give the known relation in the equidistant 
case where dj = dj + i = do and rij = n^+i = n. Unfortunately, we were not able to find 
solutions with finite values of dj ^ dj + i for both conditions. Inspired from ([6]) one may 
try 

arccos (cos(djk) — (sin(djk)). (18) 



n 



to find that this approach satisfies ffTTl) . but not fflBI) . However, choosing ( as 
in ((6]), writing dj = do(l + Aj/do), and taking the continuous limit iV — > oo with 
Aj/d -)■ u'{x) P) alters dUD to 



y 1 + u'(x) 

satisfying both ( fl6|) and (fT7|) . With the given inhomogeneous refractive index we can 
compute the electric field inside a strained, one dimensional dielectric by solving 

(dl + n 2 (x)k 2 )E(x) =0 (20) 

numerically. A comparison with the field computed via the transfer matrix method 
described in (j3J) shows excellent agreement, for sufficiently large N. 

Another way to approximate the optical field is to expand the transfer matrices 
in (jl]) for small local deformations Aj and then perform the continuous limit. This 
analytical approximation works sufficiently well for the typically small strain u' obtained 
in the scope of parameters used in this work. The lengthy results of this approach are 
presented in | Appendix B[ equation flB.6j) . 

Inserting the relation between strain and density modifications ( flOl) we finally 
obtain 

n 2 = + 1)P ~ KX) ^n 2 + (n 2 - 1) ™ ~ P , (21) 
2p - p{x) p 

where we assumed (p — p)/p 1 for the final expansion. 



2.1. Computing the reflection and transmission amplitudes 

To find solutions for the fields inside the medium with refractive index distribution 
n(x), one needs to specify initial values. As discussed for the discrete system in flH]), 
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the medium can be described in terms of a transfer matrix such that B = r\A + tD 
and Ci = tA + r r D^, if the electric fields outside the medium are given as E(x) = 
A exp(ikx)+B exp(— ikx) for x < and E(x) = Cl exp(ik(x— L))+Dl exp(— ik(x— L)) 
for x > L. The amplitudes A and Dl are determined by the intensities J; r and phases 
<f)i^ r of the fields incident from the left and right and the displacement u, as 

A = JK e ^ku(0) d D = e i<h -ik(L+u(L))_ . 22 ) 

u y £ o c V e ° c 

Therefore, the initial conditions for solutions of ( )20l) are E(0) = A + B and E'(0) = 
ik{A -B ). 

But obviously, the reflection and transmission coefficients n, r r and t strongly 
depend on the refractive index n(x). To calculate those one can either use some 
approximations, cf. Appendix B , equation (1B.5|) . or solve the field equation (1201) for 



specially chosen boundary values, e. g. 

£W(0) = iEW(0)/k = W L =► E®(L) = D L {r v + 1), E^'(L) = \kD L {r v - 1), 

E^(L) = -iE^'(L)/k = tA E^(0)=A (l + n ),E^'(0) = ikA (l + n), 

allowing the easy computation of r\, r T and t. 

It is easy to see that if n(x) is symmetric, i.e. n(x) = n(L — x), x G [0, L], then 
a beam entering from the left experiences the same medium as one from the right and 
hence r\ = r r . Note that for the homogeneous case where u' = and n = n, we recover 
the usual [XT] 

t ~ 2n 
h 2n cos(nkL) — \{n 2 + 1) sm(nkL) ' 

i{n 2 — 1) sm(nkL) 

2ncos(nkL) — i(n 2 + 1) sm(nkL) 

3. Light forces in an inhomogeneous medium 

In general, the total electromagnetic force on an object embedded in vacuum is given 
by M 

F a =/E T ^VA ( 25 ) 

J a p 

where A denotes the surface of the object, n is the normal to A and T a /3 is the Maxwell 
stress tensor 

T a/ g = £ EaE/3 + J-Bq.B/3 — ^5 a ^(e E 2 + ^-B 2 ). (26) 

Using two planes orthogonal to the direction of propagation (i. e. the x-axis) as 
integration surfaces and the plane wave fields defined in ([2]), the time-averaged optical 
force per area (pressure) on the j th slice simply reads [19] 

^• = f(i^r+i^r-i^i 2 -i^i 2 )- (27) 
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Figure 3. Illustration of the influence of a given strain u' (green) on the optical force 
density (blue) and the local intensity (red). The dotted lines show the homogeneous 
case u' = 0, the continuous lines represent the perturbations computed using the 
refractive index distribution n(x) (purple). The reflection and transmission amplitudes 
here change from r h ~ -0.29 - 0.29i, t h ~ -0.64 + 0.59i to n ~ -0.20 - 0.21i, 
r r — —0.15 — 0.24i, t ~ —0.76 + 0.51i. Please note that in this case the strain is not 
chosen such that it balances the optical forces as in (f39|) . 
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Following the beam splitter relation in ([3]) we rewrite Cj = (1 + K)Aj + K^j an d 
Dj = -i(Aj + (1 - i()Bj to obtain 



(28) 



Taking the naive limit limjv^oo Fj would give a vanishing force per slice as liniTv^oo C = 0, 
cf. (jSJ) with do = L/(N — 1). But assigning each slice to one iV th of the object's total 
length L we can define a force density J-{x) := liniAr^oo NFj/L and use 

lim = k^ 1 . (29) 

Following the derivation of the inhomogeneous refractive index (fT9~|) we can replace 
Aj + Bj = Ej^Xj+i) = E nj {\ Xj + i) — > E{x) (ITS]) and in the continuous limit it is 
reasonable to set Aj — Bj = —iE'j{^Xj + i)/k — > —iE'(x)/k, if E(x) is a solution of the 
wave equation (|2TJ|) . Hence we obtain the local optical force density 

JP( X ) = fR e [(n 2 - l)E(x)(E'(x))*\, (30) 

where we used the algebraic limit theorem for lirnjv-Kx) ^V"|C| 2 /-^ = 0- Again, the results 
from the formula above are fully consistent with forces computed using a large but 
finite number of slices (jl]) and (l2Tj) as well as an analytical approximation presented 
Appendix B, equation (IB. 9ft . 



in 



In figure [3] we compare the intensity and optical forces in a medium with 
homogeneous refractive index n to fields and forces obtained by solving f[2"Uj) and (I3"0"|) . 
respectively, for a given strain v! . 
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3.1. Identification of radiation pressure and dipole force components inside a 
homogeneous dielectric 

In the case of a medium with uniform refractive index n, i.e. with u' = 0, the force 
density computed from ( 130]) can be identified with established expressions for optical 
forces on dielectric test particles. There, the time-averaged force on a dipole at position 
x in the external field E(x,t) = Re{E(x) exp(— iut)} with E(x) = \E(x) \ exp(— iip(x)) 
reads [20] 

F L = \d x \E{x )\ 2 Rea - ±|£(x )| 2 <9^(x )Ima , (31) 

where a is the polarizability of the dipole. The first term proportional to Rea is often 
referred to as dipole or gradient force, the (dissipative) term proportional to Ima is 
called radiation pressure or scattering force [20] . 

Inside a homogeneous dielectric we may write the spatial component of the electric 
field as E(x) = G exp(ink(x — xo)) + H exp(—ink(x — xo)), where the amplitudes G 
and H are chosen such that Fresnel's conditions at the boundary of the object are met. 
Using this field and rewriting (13T]) into a force density on particles with volume density 
r/y located at x leads to 

F L (x ) = -r] Y klm{GH*}Re{an*} + §r/ v (|Gf - \H\ 2 )lm{an*}. (32) 

As the coupling parameter in ([1]) is defined as ( = kr]Aa/(2eo), with t]a denoting the 
areal particle density in each of the iV slices, we may write using fl6]) and (|29|) 

?7va = lim t^a" = lim 2-^e C = e (n 2 — 1). (33) 

This relation also resembles the Lorentz-Lorenz relation for the case of a thin gas [T7] , 
where the individual dipoles do not directly interact with each other. 

It can easily be checked that the force density in ( 1321) together with the Lorentz- 
Lorenz relation (133]) gives exactly the same result as the force computed from ( 130]) . 
if we insert the same field for a homogeneous dielectric. This demonstrates that our 
approach to calculate fields and forces from multiple scatterers is consistent with well 
known results derived from more general assumptions. 

3.2. Integrated force and trap formation 

To compute the total force on an extended dielectric in a standing wave we can use the 
same derivation as for the force on an infinitesimal slice in ( 127]) and get 

F tot = ^ (\A \ 2 + \nA + tD L \ 2 - \r r D L + tA \ 2 - \D L \ 2 ^ , (34) 

with A and Dl denoting the amplitudes at the object's left and right boundaries, as 
given in ( |22l . Defining the position of the centre of the object xo = (^(0) + L + u(L))/2, 
we can express the total force in terms of £ := xq + (0i — <j) T )/(2k) 

Ftot(0 = ~(hSl - IrSr + 2^^/7^008(2^ + V)), (35) 

with si = 1 + |n| 2 - \t\ 2 , s r = 1 + |r r | 2 - \t\ 2 , v = \r { t* - r*t\, ip = arg(nt* - r*t). 
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If Alilj-v 2 > (I\S\ — I T s r ) 2 , the force vanishes at every position £ G S + U 2_ where 
S + :={?f +fo,mGz} and S_ := { ^± - m G Z 



£o := ^( arccos 



(36) 



2«viiir 

Using some trigonometric properties one can show that stable trapping positions are 
those defined in the set H + . For £ G 2_ we find i^ot^o) > and hence S_ is a set of 
unstable trapping position. Linearising the total force in ( 13 5 p around stable trapping 
positions £ G ^+ leads to a trap stiffness of 

k = f v^i^ 2 - (/ lSl - / rSr )2. (37) 

Since the reflection and transmission coefficients strongly depend on the object's size, 
one finds that the parameter if> can change its sign abruptly for certain values of L. This 
leads to sudden jumps between a low- and high-field-seeking behaviour for a trapped 
object [2U E2} E3j E]. Figure [5] shows how trap position and trap stiffness is changed 
by the strain induced on the object by optical forces. 



4. Self-consistent balancing of optical force and elastic back-action 

In the previous chapters [2] and [3] we found expressions for the local fields and the 
local optical force densities in deformed, dielectric media. But depending on the given 
elastic properties, the strain will result in stress which typically tries to compensate the 
external volumetric forces. 

In this chapter we will investigate the behaviour of a linear elastic, dielectric object 
subjected to the optical forces described by ( 1301) . More precisely, we will provide a 
framework to compute the equilibrium configuration between the optical forces and the 
elastic counter reaction in a self-consistent manner. In our computations we will assume 
only optical forces and neglect thermal or piezoelectric effects as well as surface tension. 

Mechanical equilibrium between some general volume force density / and the 
resulting stress denoted by the tensor a is given by Cauchy's equilibrium equation [16] 

/i + 5>i^ = 0, (38) 

3 

for i,j denoting the coordinates of the system. Since the model discussed here considers 
only one relevant dimension, this equilibrium equation simplifies to / + d x a = 0. The 
constitutive relation for a linear elastic, one dimensional object simply reads a = £u', 
where S is Young's modulus and u' is the local strain [T6] . 

Hence we see that an equilibrium between the optical force density and the elastic 
strain requires 

T{x) + £u"{x) = 0, (39) 

at every position x G [0,L], with T being a solution of ( 130]) . Note that the electric 
field computed from ( 120|) also depends on the amplitudes at the edges of the object and 
therefore also on the displacement u, cf. ( )22l) . 
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Solving (j3HD for an equilibrium requires boundary conditions on the displacement 
u and the strain u' which are determined by the given setup. The displacement is fixed 
by the assumed trapping mechanism, e. g. if the object is trapped by a standing wave, 
we have to fulfil (u(0) + L + u(L))/2 G E + . But note that £ depends on the reflection 
and transmission coefficients and hence also on the deformation u', cf. (1361) . 

The strain has to be chosen such that the stress a = Su' at each surface balances 
external surface forces [16J. Assuming for the moment an object subjected to volumetric 
optical forces only, we integrate the equilibrium equation ( 1391) at obtain 

0= f (F(x) + o-'(x))dx = F tot + a(L)-o-(0). (40) 
Jo 

For an object trapped by light fields, we get F tot = and er(0) = a(L) = 0, due to 
the lack of surface pressure. In chapter 15.21 however, we fix the slab by an external 
mechanism balancing the total optical force via surface interaction. Hence if the left 
boundary of the slab is retained at x = (i. e. u(0) = 0), then o~(L) = and cr(0) = F tot . 

To solve equation (1391) numerically, we use an iterative approach where the 
equilibrium condition is rewritten in the form 

F{x)[u h u[\+£< +l {x) = ^ (41) 

with Ui and u[ denoting the displacement and strain obtained by the i th iteration step 
and T(x) [u^ u^] is the force density computed using Uj and u\. The updated u' i+1 and 
u,i + i can then be obtained by simple numerical integration, with integration constants 
chosen in accordance to the boundary conditions of the used setup. With the updated 
optical force densities J-(x)[ui + i, u' i+1 ] one can compute the next step of the iteration. An 
obvious choice for initial values is a homogeneously shifted distribution (i. e. a constant 
u Q ) with a given starting length L and refractive index n. 

This iterative scheme proved to be sufficiently exact but significantly faster 
compared to other methods of solving nonlinear equations, such as Newton's method. 
We also confirmed our computations with force densities obtained from the transfer 



matrix approach in and the analytic approximation described in |Appendix B 
respectively. 



5. Examples and physical interpretation 

In our basic considerations above we always assumed the object to be exposed to two 
counterpropagating laser beams of the same, linear polarisation forming a standing wave. 
These results can easily be extended to describe situations with only one incident beam 
or with two counterpropagating beams of different polarisation. In the latter case, one 
has to calculate the intensities and forces separately for each polarisation direction. 

In this section we will present four showcase examples to give insight into the large 
variety of possible results. The first two examples deal with the case where an object 
is trapped by two counterpropagating beams. The latter examples will treat the case 
where the object is illuminated by only one beam and externally fixed at one end. For 
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Figure 4. Relative length change, transmission \t\ 2 and reflection \r\\ 2 for a dielectric 
slab trapped by two plane waves forming a standing wave (blue) or having orthogonal 
directions of polarisation (red lines). The incoming intensities and Young's modulus 
are related as I\ = I x = 0.05£c, where c is the speed of light in vacuum; the refractive 
index is chosen as n = 1.3 + 0.0025i. The discontinuities in the standing-wave case 
stem from abrupt jumps in the stable trapping positions at L — mA/(2Re[n]), to e N 
(grid lines) . The circles indicate the values used for the examples in figure [6] and [7] 



each of the given examples one has to specify the boundary conditions for u and u', as 
discussed in chapter 0J 

For all considered setups we will see that the interaction between optical forces and 
elastic back-action strongly depends on the ratio between the initial length L and the 
wavelength of the deformation beam in the unperturbed medium, A/Re[n]. Concerning 
the intensity and the elastic properties of the medium we find that all results grow 
linearly in (ij + I r )/(Sc), at least in the scope of parameters where a solution of the 
equilibrium equation ( |39|) could be obtained within a reasonable error tolerance. That 
is why the local intensities, J(x) = Eqc\E(x)\ 2 /2, and force densities are given in units 
proportional to £ in the upcoming figures. Note that the numbers used in the simulations 
are unrealistic in order to exaggerate the effects, since an intensity of I = O.ISc would 
imply I ~ 30W//im 2 for £ ~ IMPa. 

5.1. Example: Object trapped by two counterpropagating beams 

As argued above, an object with initial length L subjected to optical forces will in general 
experience local deformations and an overall length change. Figure H] shows the relative 
length change AL = (L — L)/L for different initial lengths L in the cases where an object 
is trapped in a standing wave (blue lines) or by two beams of orthogonal polarisation (red 
curves) and equal intensity. Obviously, both configurations are symmetric regarding an 
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Figure 5. Trap stiffness k (|37|) and stable trapping position £o € (|36[) for an 
object trapped in a standing wave. The same parameters are used as in fig. 21 
i.e. I\ = J r = 0.05£c and n = 1.3 + 0.0025L A comparison of the blue curves 
denoting equilibrium solutions of (|39p with the orange, dashed lines computed for the 
homogeneous objects shows that the deformation induced by optical forces significantly 
increases the total trapping strength. For the given case where I\ = 7 r , the trap 
positions remain unchanged. 



inversion of x at the centre of the object and therefore we find r\ = r r , as mentioned in 
section 12.11 

Surprisingly, in the standing wave case we observe abrupt switching from strong 
compressive to stretching behaviour around certain initial lengths. A comparison with 
chapter [3.21 and earlier discussions in [TT] shows that these switches are concurrent with 
jumps of the stable trapping position £ - Generally speaking, objects with small values 
of kL are trapped at local maxima of the intensity in the standing wave. But for larger 
objects, the term if) in (|36|) abruptly changes its sign and the object seeks positions 
centred around intensity minima. As indicated by the dotted grid lines, these jumps 
occur at lengths of minimal reflection |rh| 2 and maximal transmission |th| 2 (EIJ), i- e. at 
L = mA/(2Re[n]), m G N. 

For larger objects we observe a slight decay of the maximal relative length change, 
which is found also for computations with Im[n] = and thus is not caused by additional 
radiation pressure only. However, a glance at the reflectivity and transmission shows 
that the self consistent deformation prevents configurations with zero reflectivity which 
would usually result in maximum elongation or compression. For the given standing 
wave trap we generally observe that deformations computed with the steady state 
equation ( 139|) tend to increase reflectivity and decrease transitivity compared to a 
homogeneous medium, for both lm[n] = or lm[n] ^ 0. 

In figure [5] we compare the trap stiffness k from fl37j) and the trapping position £ 
from ( |36l) between unperturbed and self-consistently deformed objects. We find that 
the deformation significantly increases the trap stiffness, even if the total size of the 
object remains unchanged, cf. figure HI An elastic object in a standing wave therefore 
assists in enforcing its own trap, just like a rabbit who starts to dig when captured in 
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x [in units of A] x [in units of A] 

Figure 6. Dielectric objects trapped in a standing wave configuration with I\ = 
I r = 0.05£c. The initial lengths are chosen as L ~ 0.19A (a) and L ~ 0.4A (b); 
n = 1.3 + 0.0025i. As mentioned in (fT0|) . positive strain u'(x) (green curves) denotes 
larger distances between the particles and hence a reduction of material density. The 
density is increased at positions where the force density (blue lines) T is zero and 
J-' < 0, i. e. where the light intensity (depicted red) has a maximum value. As indicated 
by the marks in figure H] and verifiable from the signs of u'(x), the object on the left 
gets squeezed whereas the right hand side example shows a stretched object. 

a pit. Note that for a compressible slab in a standing wave there is no critical length of 
zero trap stiffness, but deformation always leads to a stable trapping position. 

Two examples for optical force densities and the associated deformation in the 
standing wave setup are shown in figure EJ There we see that the strain is negative 
(i. e. the material density is increased) at positions where the force density changes from 
positive values, denoting a force pushing to the right, to negative values associated with 
a force pushing to the left. One can clearly see the difference between the compressive 
situation (a) where the object is trapped at maximal intensity and situation (b), where 
the trapping occurs at minimal intensity and the object experiences a stretching force. 

To trap a dielectric slab with two non-interfering plane waves of orthogonal 
polarisation, the intensities of said beams have to be equal, i. e. I\ = I r . In this case there 
are certain starting lengths L ~ (2m + l)/(4Re[n]), m G No, for which the intensities 
inside the unperturbed object add up to a constant value. Hence for these specific 
lengths the dominant gradient forces add up to zero and the object's length remains 
unchanged. Apart from that we find an expanding behaviour for L > A/(4Re[n]), as 
can also be seen in the examples in figure [7J 

In the case of non-interfering beams, figure H] also shows that the reflectivity and 
transitivity is no longer periodic in the object's length L. For the reflectivity we find 
that the zeros are shifted from L = mA/(2Re[n]), m G N, towards smaller lengths and 
no longer coincide with lengths of maximal elongation. 
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Figure 7. Dielectric objects trapped by two counterpropagating beams with 
orthogonal polarisation and I\ = I T = 0.05£c. As in figure [6j the initial lengths 
are L ~ 0.19A (a) and L ~ 0.4A (b); n = 1.3 + 0.0025L As can be seen from figure |H 
the length of the first example is chosen such that the intensities (red curves) inside the 
medium almost add up to a constant value. Hence the force densities (blue curves) and 
the self-consistent deformation (green lines) practically vanish (please note the changed 
scale on the axes). For the second example we obtain a stretching of AL ~ 6.6 x 1CP 3 . 

5.2. Example: Object fixed at the left boundary and illuminated by one beam. 

Figure [8] shows the relative length change in the case where the left edge of a dielectric 
slab is fixed by some external mechanism. Here we observe a striking difference whether 
the object is illuminated from the left (blue) or right hand side (red curves). In the 
first case we find only stretching behaviour with minor oscillations of the relative length 
change AL. However, if the object is illuminated from the right hand side (i. e. the 
beam is incident on the free surface) we again find both compression and elongation, 
depending on the initial length L. 

Furthermore we observe a stronger impact of absorption than for the previous case 
with two counterpropagating beams. Now we see that if Im[n] 7^ 0, then AL increases 
(light incident from the left, radiation pressure pushing the object to the right) or 
decreases (light incident from the right) for larger initial lengths L. As we see from the 
dash-dotted lines, the length change continues to oscillate around a constant value also 
for large values of L, if Im[n] = 0. 

For the reflectivity [transitivity] , the self consistent strain again results in a shift of 
the minimal [maximal] values to the left of L = m\/ (2Re[n]), m G N. Note that for the 
given, asymmetric setup n 7^ r r , but the difference in the total reflectivity |ri| 2 — |r r | 2 
remains very low for the given parameters. 

To explain the different results for left and right incident beams, we take a closer 
look at the interface where a single beam exits the slab, e. g. at hand of figure [9j Let us 
assume a beam entering from the left and J r = 0. Then the intensity is not only constant 
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Figure 8. Relative length change, transmission \t\ 2 and left reflection |ri| 2 for a 
dielectric slab fixed at the left edge and illuminated by only one beam incident from 
the left (blue) or right (red curves) hand side, respectively. The initial refractive index 
is chosen as n — 1.3 + 0.0025i, the dash-dotted lines show the results for the non- 
absorptive case n = 1.3 (in the third plot for \r\\ 2 the additional lines were omitted 
to avoid ambiguity). The beams are of intensity ii, r = O.lfc, the circles indicate the 
values used in figure [9] The grid lines mark the locations of the minima [maxima] of 
the reflectivity [transitivity] in the homogeneous case at L = mA/(2Re[n]), m E N. 




-0.25 0.25 0.5 0.75 1 -0.25 0.25 0.5 0.75 1 

x [in units of A] x [in units of A] 

Figure 9. Two dielectric objects fixed at the left edge and illuminated from the left (a) 
or from the right hand side only (b), with 7i ir = O.lfc. As predicted from (j42|) . the 
intensity (red) is locally maximal and the force densities (blue lines) are (almost) zero at 
the interface where light exits the slab. For left incidence we observe only positive (i. e. 
stretching) strain u' (green lines) whereas both elongating and compressive deformation 
is feasible for right illumination. 
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on the right of the object, but Fresnel's formulae also tell us (for a homogeneous object) 

£/(L) = and £l(L) = k 2 I(L)(l-2Re[n] 2 + 2lm[n} 2 )<0, (42) 

and hence the intensity has a local maximum at the edge where the beam exits the slab. 
For the usually dominant gradient force J-'gr(x) ~ 4-I{x) we thus find J-'gr(L) = and 
J-g r > left of the surface. Considering the steady state equation (1391) and omitting the 
scattering force shows that 

F{L) + Eu"{L) ~ JF ff (L) + £u"{L) = (43) 

and hence the strain v! has a local minimum at the right edge x = L. If light enters 
only from the right hand side we find analogous behaviour for the left edge x = 0. Since 
minimal strain corresponds to a maximum in local material density (11 01) . we conclude 
that the gradient force tends to accumulate material at the surface where a single light 
beam exits the object. A similar statement holds for the aforementioned case with two 
counterpropagating beams oscillating in orthogonal polarisations, but then one has to 
add up the forces generated by the two beams. 

In total, the constraints on the strain as derived from f l42|) . (1431) and chapter [4] are 
found as 

left incidence only, I T = : u'(0) = J-tot/£ > 0; u'(L) = 0, is a minimum, (44) 
right incidence only, I\ = : u'(0) = J-" t ot/£ < 0, isaminimum; u'(L) = 0. (45) 

So if I\ = 0, the strain is fixed at a minimum with negative value on the left edge, at zero 
on the right edge, and oscillates proportional to the intensity in between. So in total we 
find both negative (compressive) and positive (expanding) deformation, depending on 
the length and refractive index of the object. For I r = we see a positive strain at the 
left edge and a minimum with u'(L) = at the right boundary. Hence the deformations 
oscillate between zero and some positive value and always lead to a total stretching 
behaviour. 



5.3. Identifying length changes by probing the reflectivity 

As one can see from figures H] and El the relative length change AL obtained for a given 
ratio of total intensity to Young's modulus can be imperceptibly small, especially if 
the original length is not chosen in an optimal relation to the vacuum wavelength of 
the trapping beam, A. But one possibility to detect minor stretching or squeezing for 
arbitrary initial lengths can be found in the use of a second, weaker laser probing the 
change in the reflectivity of the medium: Assuming a non-dispersive medium, a probe 
laser with a vacuum wavelength matching the Fabry-Perot condition 

2LRen . . 

A p ~ , m G N, 46 

m 

will travel through a slab of length L without being reflected, i. e. |^h(A p )| 2 — 0. Turning 
on a powerful laser will deform the object and hence also change the reflectivity for the 
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Figure 10. Relative change in the reflectivity of a strained object for different 
probe laser wavelengths A p , for slabs trapped in a standing wave pattern with 
h = I T = 0.0005£c (red) and I\ = I T = 0.005£c (blue curves), n = 1.3 + 0.0025i. In 
the left frame we used an initial length L = 0.7 A and obtained compressive behaviour 
with AL ~ -0.2 x 10~ 3 (red) and AL ~ -2.0 x 10~ 3 (blue). The results on the 
right indicate stretching with AL ~ 0.3 x 10~ 3 (red) and AL ~ 3.1 x 10~ 3 (blue) for 
L = 0.8A. 



weak probe beam. Figure [10] shows the relative change of reflectivity 
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between the unperturbed and the strained medium for different wavelengths of the 
probe laser beam. We can see that as A p crosses values defined in (JIB"]) . 6~R changes from 
negative to positive values if the object is compressed (i. e. AL < 0) and vice versa if 
AL > 0. 

6. Estimating the deformation by computing the photon momentum 
transfer on a surface 

There exist numerous experimental and theoretical papers reporting optical stretching of 
deformable objects, such as biological cells [2HI25], or light induced outward bending of 
liquid-gas surfaces [261 127]. m the mentioned publications, the light-induced deformation 
is estimated by considering an effective photon momentum change at the transition 
from one medium to another. In this context, the optical forces emerge as surface 
forces, acting on the interface between two regions of different refractive index. Since 
the considered materials are incompressible, the refractive index in each region remains 
constant. 

Let us try here a similar approach to estimate the deformation of an elastic object 
and put our results in context with these earlier works. But note that the very different 
physical properties of a linear elastic medium as compared to incompressible water, 
plane waves instead of Gaussian beams and a wave description instead of geometric 
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optics, do not allow a straightforward comparison of the results. But nevertheless we 
can investigate whether our findings based on a volumetric description of optical forces 
are compatible with a concept of surface forces due to photon momentum exchange. 

Following the line of [2U [25j [26] one estimates the time averaged force per area on 
an interface separating two regions with different indices of refraction n\ 7^ n 2 G M and 
fields E(x) = Aexp(inikx) + B exp(— m\kx), x < 0, and E(x) = C exp(in 2 kx), x > 0, 
as 

Ap l mc ( Pl (l + R)-p 2 T) 

Fl ' 2 = At = Wc ' (48) 

Here J inc = nice \A\ 2 /2 is the total energy-flux density entering the system, R = 
\B\ 2 /\A\ 2 and T = n2|C| 2 /(ni|A| 2 ) is the reflected and transmitted fraction the energy- 
flux and pi = Hkrii describes the momentum of a single photon in a medium with index 
Hi, i = 1,2. As in [2H [221 EE] we here used Minkowski's version of the momentum 
of light in dielectric media. Out of curiosity about simple but puzzling arguments 
on stretching or compression of media in connection with the Abraham-Minkowski 
controversy [12], [26l [281 I2H] we also included results computed by naively inserting 
Abraham's result for the photon momentum, pi = hk/ni, in figure [TTJ One must note, 
however, that Abraham's stress tensor would also require a material tensor component. 
A thorough calculation should always give the same results, independent of the used 
version of stress tensor [28, 29J. 

Neglecting internal reflections, the force on an extended object with index n 2 
embedded in a medium ri\ and subjected to a single beam is estimated to give 
Ftot = Fi j 2+TF2 t i- The deformation of such an object is then simply the difference in the 
two forces on the surfaces, reading F def = TF 2 i — Assuming a linear elastic medium 
with Young's modulus S, the relative length-change can be estimated by AL = 
For the values used in the previous examples n\ = 1, n 2 = 1.3, Iq = 0.1£c, we obtain 
AL ~ 0.059 when using Minkowski's momentum and AL ~ —0.046 for Abraham. 
These deformations have about the same order of magnitude as our full self-consistent 
computations, depicted e. g. in figure E] but do not depend on the length of the object. 

Formally one can refine these calculations from ( HHl) and include also light incident 
from the right such that E(x) = C exp(in 2 /ca;) + D exp(— in 2 kx), x > to obtain 

Fl ' 2 = W* 1 - 4 ' 2 + |B|2) - W^ 12 + |c|2) ' (49) 



For D = this reduces to ( 1481) and for n\ = n 2 = 1 and pi = p 2 = hk we recover the 
force derived previously with the Maxwell stress tensor f [27|) . This allows to formulate 
a generic wave optics extension for the deformation estimated above in the scope of 
geometric optics, now including also size dependent reflection and transmission. The 
total deformation pressure on an object with length L and homogeneous refractive index 
n G 1R surrounded by vacuum then reads 

Fdef = ^(iGf + |# I 2 ) - |(|A| 2 + \B\ 2 - |Cf - \D\ 2 ) (50) 

where the amplitudes inside the medium are computed using Fresnel's relations G = 
{{n + I) A + (n - l)B)/(2n), H = ((n - I) A + (n + l)S)/(2n), the amplitudes 
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Figure 11. Expected deformation due to the change in photon momentum at the 
interfaces between vacuum and a homogeneous dielectric of refractive index n = 1.3, 
as obtained from (|50p. For the blue lines we used Minkowski's photon momentum 
p = hkn, the red curves are computed using Abraham's p = hk/n. The left figure 
shows the situation where a homogeneous object is trapped by a sanding wave with 
I\ = I x = 0.05£c. As in figure IH the deformation changes abruptly, if the object 
switches from a low-field to a high-field seeking behaviour. The right hand figure 
depicts the relative length change for an object illuminated from only one side with 
I\ = 0.1£c, / r = 0. In both situations we find that Minkowski's momentum only leads 
to stretching, whereas the naive adaption of Abraham's photon momentum would 
solely result in compression. 



outside are connected by the homogeneous reflection and transmission coefficients ( 1241) . 
B = r^A + t^D, C = t^A + r^D, and the incoming A and D are given in fl22l) . As 
expected, a similar calculation for the total force F to t gives the same result as we obtained 
previously in (I34I) . 

The resulting relative length change AL = F^ c i/£ is presented in figure [HI There 
we find that the estimations using optical surface forces even qualitatively differ from 
the results we obtained with the present description using the full, volumetric optical 
forces, cf. figures H] or El even if the force on the surface is adapted to include interference 
due to internal reflections. 

An intuitive example is the object of length L = X/(2n) where r^ = and 
th = 1- In a standing wave trap with /; = I r this object is then trapped at 
£o = + (0i — 4>r)/(2k) = 0, cf. chapter [3721 or figure [U Using Minkowski's p n = hkn 
we then find that the force due to photon momentum transfer vanishes at each surface 
and hence AL = in figure [TTJ 

But from the examples in figure [6] we deduce that the intensity at the surface is zero, 
yet the object will contract due to the intensity maximum at its central position. This 
is because the dipole force pulls each volume element towards the next local maximum 
of intensity, regardless of whether this volume element is located at the surface or in the 
bulk of the medium. 
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In fact, none of our calculations or simulations showed any distinctive effects 
suggesting a surface force at the boundaries of a dielectric [IT]. This is supported 
by computations on large but finite stacks of polarizable slices where the forces on 
the first or last slice qualitatively do not differ from those on the second or next 
to last, respectively We therefore conclude that optical forces have to be treated 
as real volumetric forces |T2l [30] and that a description using the change of photon 
momentum at the surface of a medium is inappropriate, regardless of using Abraham's 
or Minkowski's momentum. 

7. Conclusions 

Using an implicit calculation of optical fields and forces allows to self-consistently 
determine the stationary local deformations of an elastic object, where the local stress 
balances the local light forces by elastic back action. These solutions show a surprisingly 
variable and nonlinear dependence on the chosen parameters. Generally we see a length 
and illumination dependent, spatially quasiperiodic strain pattern, which can lead to 
length stretching as well as compression. As expected, standing wave configurations 
yield the strongest forces and effective length changes with a clear resonant structure 
for special ratios of initial object length L and trap beam wavelength A. In the standing 
wave setup, variations in the trap wavelength A lead to discrete jumps of the stable 
trapping positions. At L — mA/(2Re[n]), m e N, the particle switches from a position 
centred around an intensity maximum to one around a field node, which is associated 
with changes from compression to elongation of the object. Interestingly, in particular 
close to these instability points, this generally leads to an increase in trap stiffness. 
We expect that this indicates possible bistability between high and low field seeking 
behaviour for certain lengths very close to L = mA/(2Re[n]). 

Although the calculations presented here were performed in the scope of elastic 
media, we believe that the model can be extended to deformable but incompressible 
media like water or even dilute gases. Here in particular stability thresholds for the 
homogeneous solutions should prove physically very interesting, as they could lead to 
stationary flows, periodic density oscillations or light induced density pattern formation 
and particle ordering in a gas. 

Here we limited our considerations to the case, where a steady state solution can 
be found. As for other nonlinear dynamical effects [H], there are of course regions in 
parameter space, were no stationary solutions exist and we find self sustained oscillations 
or even disintegration of the material. Indications of this behaviour appear e. g. in a 
non converging iteration procedure. At this point we leave this to future work. 
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Appendix A. General features of transfer matrices 

In equation (jHJ) we already used the concept of a transfer matrix to couple the plane 
wave amplitudes left and right of a beam splitter. Let us generally define the set of 
transfer matrices as 

T := {t G C 2X2 | 3r h r 2 ,t eC:T = j(* T A = T(r 1; r 2 ,t)}. (A.l) 

One can easily show that Tj • T 2 G T for all Ti, T 2 G T and • here denoting the usually 
omitted matrix multiplication. Since also 

[T(r u r 2 ,t))-^-( 1 ri )gT (A.2) 

t V _r 2 * - r l r 2 J 

we conclude that (T, •) is a group. 

To motivate definition (lA.ljl . let us consider two plane waves E\(x) = Aexp(ik(x)) + 
B exp(— ik(x)), x < 0, and E T (x) = C exp(ifc(x — L)) + D exp(— ik(x — L)), x > L, left 
and right of a dielectric with a homogeneous refractive index n. Hence their amplitudes 
are coupled as 

(A.3) 



where S„ ljn2 couples the amplitudes at the intersection from a region with refractive 
index n\ to a region with index n 2 and denotes the propagation matrix over a 
distance d, i. e. 

1 (n 2 + n l n 2 -n l \ ( e ikd \ 

s ni,n 2 = 7^— , and P d = . fcd . (A.4) 

2n 2 \n 2 - n x n 2 + n x J \ e XKa J 

One can easily show that S n> i ■ P„l ■ Si jTl = T(r"h, Th, th)> with rh and th as given 
in ( 124"1) . But the attempt to construct a total transfer matrix for stacked media such as 
Sn K ,i-Pn K d K -S nK _ ltnK ■ • • S n2>ril • P nidl • S lm with arbitrary n*, d h i = 1, . . . , K, will show 
that one reflection coefficient is not enough. However, such a system can be described 
by a more general T(r 1; r 2 , t) G T such that S = + tD and C = r 2 .D + tA. 

Appendix B. Analytical approximations for electric fields and forces for 
small deformations 

In section [2] and in an earlier work [11] we showed that for equally spaced slices 
the amplitudes of the electric fields are related as (Aj+i, Bj + i) T = TV(Al, B\) t , 
with T h := P df) M B s- Choosing the coupling ( as in (jSJ), the eigenvalues of T h read 
exp(±mfcc/o). This leads to 

lim T h ^ =: T(x) = - ( + " 9 f{ ) (B.l) 

at->oo yj 2n\f(-x)-g(-x) f(x)+g(x)J 1 ; 

where x = linijv^oo L(j — 1)/(N — 1) and 

f(x) = ncos(nkx) — isin(nkx), ^ 

g(x) = n cos(nkx) — in 2 sm(nkx). 
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Allowing small local density variations dj — d = Aj we may use ~ P^ + ikAjO~ z Pd 
to expand the relation between the field amplitudes from (j3J) as 

( B +1 ) ~ [ ThJ + ik ^ A - T h J " m ^ T h m ] ( g 1 ) , (B.3) 

\ i+l / m=l ^ ' 

where = diag(l, —1). In the limit of infinitely many slices within a finite length L, 
the sum above can be rewritten to an integral and we obtain 

( B{x) ) " [ T(X) + ik [ U ' iy)T{X ~ y)azT{y)dy ] ( B% ) • (R4) 
Since limjy^oo M B s = id2, the total reflection and transmission coefficients for the 
displaced object can be read off and expanded in linear order of u' from equations (J5]) 
and flRlj) 

-L 



* - + wD^wy / u '(y)(f(y)f( L -v) + a{y)g{L - y))d y , 

J 

n * r h + (/(L) % (£))a fn'{y){f{L-yf-g{L-yf)dy, (B.5) 

*/ 



rh + TTri^nw *'(v)(/(v) -0(i/) a )dy- 



(/(L)+g(L))2 

Here and th denote the reflection and transmission amplitudes of a homogeneous 
medium with length L and refractive index n, as given in f )24|) . As mentioned in 
section 12.11 we find that symmetric local strain, i.e. u'(x) = u'(L — x), x G [0,L], 
results in n ~ r r and antisymmetric strain gives t — th and r\ — r^ ~ rh — r r . 

Using the reflection and transmission amplitudes from (IB. 51) to replace -B(O) = 
r\A(0) + tD(L) in (IB .41) leads to analytical approximations for the local field amplitudes 
inside a deformed medium 

u'(y)a x (x,y)dy + / u'(y)a h (x,y)dy, 

J \ (B.6) 

rx i>L v ' 

B(x)~b h (x)+ u'(y)b x (x,y)dy+ u'(y)b L (x,y)dy, 
Jo Jo 
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with the individual terms reading 

A (f(L -x)+ g(L - x)) + D L (f(x) - g(x)) 



a h {x) 
b h (x) 
a x (x,y) 

K(x,y) : 

ah(x,y) 



f(L) + g(L) 

MfjL -x)- g(L - x)) + D L (f(x) + g(x)) 
f(L)+g(L) 

A (f(L - y)f(y -x) + g{L - y)g(y - x)) + 
+D L (f(y)f(y -x)- g{y)g{y - x))j , 
A (f(L - y)f(x -y)- g(L - y)g(x - y)) + 
+D L (f(y)f(x -y)+ g{y)g{x - y)) 

ik(f(x) -g(x)) 
n{f{L)+g{L)f 

ik(f(x)+g(x)) 



ik 

n(f(L)+g(L)) 



-ik 

n(f(L)+g(L)) 



(B.7) 



{A,U(L-yf-g(L-yf) + 

+D L {f(y)f{L-y)+g(y)g{L-y)) 
^A (f(L-y) 2 -g(L-y) 2 ) + 



n(f(L)+g(L)Y 

+D L (f(y)f(L-y)+g(y)g(L-y))^ 

The functions / and g are given in ( IB. 21) . the amplitudes at the boundary A = A(0) 
and Dl = D(L) are determined by the incoming intensities and the displacement and 
can be read off from (|22l) . 

To obtain an approximation for the forces we use (}28l) . take the limit ^{x) : = 
limjv^oo NFj/L with lim A r_ 5 . 00 N(/L = k{n 2 — l)/2 and insert the amplitudes from (IB. 61) 
to find 

JF( X ) = ^ilm{(n 2 - l){A{x) + B(x)){A(x) - B(x))*} 

px pL 

J c (x)~T h (x)+ u'(y)J r x (x,y)dy+ u'(y)J r L (x,y)dy, (B.8) 
Jo Jo 



where 
T h (x) 

T x (x) 



|(/(L 2 ) ^' (L)| J m{(n 2 - l)(A f(L - x) + D L f(x))(A g(L - x) - D L g(x))*} 
^Im{(n 2 - l)[(a h (a;) - b h {x))*{a x {x,y) + b x {x,y)) + 
{a h (x) + b h (x))(a x (x,y) - b x (x,y))* 
Fl(x) = *flm\(n 2 - l)\(a h (x) -b h (x))*(a L (x,y) + b L (x,y)) + 



(B.9) 



(a h (x) + b h (x))(a L (x,y) - b L (x,y)Y 

with the terms ax and &x as given above in (1B.7|) . 

Just like the field amplitudes in ( IB. 61) and the reflection and transmission coefficients 
in (1B.5j) . the above result is only an approximation for the case of small deformation 
u' . That is, we neglect products of A;A m in (IB. 31) and hence also correlations of 
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type J J u'{yi)u'{y2) ■ ■ . dyidy 2 - But a comparison of the approximated results in the 
continuous limit with solutions of the wave equation (120]) or numerical computations for 
a large but finite number of beam splitters confirmed that this first order expansion is 
sufficient for the scope of parameters used in this work. 
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